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Using high statistics samples of charged current interactions, MiniBooNE reports a measure- 
ment of the single charged pion production to quasi-elastic cross section ratio on mineral oil (CH2), 
both with and without corrections for hadron re-interactions in the target nucleus. The result is 
provided as a function of neutrino energy in the range 0.4 GeV < < 2.4 GeV with 11% precision 
in the region of highest statistics. The results are consistent with previous measurements and the 
prediction from historical neutrino calculations. 



Future neutrino oscillation experiments will operate in 
the 1 GeV energy range, where charged current quasi- 
elastic scattering (CCQE, i/^n Ii'~p) and charged cur- 
rent single pion production (CCl7r+, V/^iX fi^ir^X') 
are the dominant interactions. Because such processes 
are the largest contributors to the event samples in such 
experiments, there has been much interest in making bet- 
ter determinations of their cross sections. At present, 
the ratio of CCl7r+/CCQE cross sections has been mea- 
sured to ~30% precision based on small event samples 



[IlEJIS]. A high statistics measurement of these processes 
necessarily requires the use of nuclear targets where final 
state interactions obscure the actual value of the ratio 
of the cross sections on nucleons. Experimentally, it is 
the cross section on complex nuclei including the effects 
of final state interactions which is more relevant (exper- 
iments can only identify particles that actually exit the 
struck nucleus). In this letter, MiniBooNE reports the 
first measurement of the observed CCl7r+/CCQE cross 
section ratio as a function of neutrino energy including 
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the effects of hadronic re-interactions. Additionally, an 
underlying ratio at the nucleon level is extracted to facil- 
itate comparison with prior measurements [TJ |3] . Precise 
knowledge of this cross section ratio is particularly im- 
portant for future disappearance searches, in which 
CCl7r+ events typically constitute either a class of sig- 
nal events or a large background to the CCQE signal. 
The uncertainty on the CCl7r+/CCQE cross section ra- 
tio therefore limits the precision of these measurements. 

The Booster Neutrino Beam at Fermilab provides a 
neutrino source which is particularly well-suited to make 
this measurement; about 40% of j/^ neutrino interac- 
tions in MiniBooNE are expected to be CCQE and 24% 
CCl7r+. The beam itself is composed of 93.6% i^^ with 
a mean energy of about 800 MeV and 5.9% (0.5%) i>p 
(ve) contamination [l]. The neutrinos are detected in 
the MiniBooNE detector a 12.2 m diameter spherical 
tank filled with 818 tons of undoped mineral oil located 
541 m downstream of the beryllium target. At the en- 
ergies relevant to this analysis, the products of the in- 
teractions produce primarily Cerenkov light with a small 
fraction of scintillation light [5]. The light is detected 
by 1280 8-inch photomultiplier tubes (PMTs) which line 
the MiniBooNE inner tank. This inner tank region is 
optically isolated from a surrounding veto region, instru- 
mented with 240 PMTs, that serves to reject incoming 
cosmic rays and partially contained neutrino interactions. 

Neutrino interactions within the detector are simulated 
with the v3 NUANCE event generator |6]. CCQE in- 
teractions on carbon are generated using the relativistic 
Fermi gas model 7 tuned to better describe the observed 
distribution of j/^ CCQE interactions in MiniBooNE [8J. 
Resonant CCItt^ events are simulated using the Rein 
and Sehgal (R-S) model as implemented in NUANCE 
with an axial mass M)^ = 1.1 GeV. The angular distri- 
bution of the decaying pions in the center of mass of 
the recoiling resonance follows the hclicity amplitudes of 
0. In MiniBooNE, 87% of CCl7r+ production is pre- 
dicted to occur via the A(1232) resonance, but 17 higher 
mass resonances and their interferences, as well as a non- 
resonant background 9 that accounts for roughly 6% of 
CCl7r+ events, are also included in the model. Coher- 
ently produced CCItt'*" events are generated using the 
R-S model [10 with the R-S absorptive factor replaced 
by nuance's pion absorption model and the overall 
cross section rescaled to reproduce MiniBooNE's recent 
measurement of neutral current coherent tt*^ production 
[TT] . Coherent 7r+ production is predicted to compose 
less than 6% of the MiniBooNE CCl7r+ sample due to 
the small coherent cross section [T2 [T^] and the domi- 
nance of the A++ resonance. A GEANT3-based detector 
model simulates the response of the detector to par- 
ticles produced in these neutrino interactions. 

To select charged current interactions, simple re- 
quirements on the amount of charge detected in the tank 
(> 175 tank PMT hits) and in the veto region (< 6 veto 



PMT hits), location of the event in the tank (< 500 cm 
from the center of the detector), and event time (event 
must occur while the beam is passing through the detec- 
tor) are first applied. Further requirements on the num- 
ber of decay electrons in the event are then used to isolate 
CCQE from CCl7r+ interactions, i^f^ CCQE events are 
selected by requiring the detection of a single electron 
(from the decay of a stopped muon) within 100 cm of 
the endpoint of the muon track |H] . Identification of the 
decay electron is possible because it follows the detection 
of the muon by a distinct time interval, CCItt"*" in- 
teractions are identified by requiring the detection of two 
electrons (from the decay of the muon e~z/^z7g) 
and pion (7r+ l^'^v^, ^+ e^v^Ve)), at least one 
of which must be within 150 cm of the endpoint of the 
muon track. The model dependence of the event selec- 
tion is rather small since we require only that the /i~ 
and 7r+ decay. After cuts and with 5.58 • 10^° protons on 
target, the CCQE data sample contains 193,709 events 
and the CCl7r+ sample 46,172, making these the largest 
samples collected in this energy range by more than an 
order of magnitude. 

The CCQE and CCl7r+ reconstruction requires a de- 
tailed model of light production and propagation in the 
tank to predict the charge distribution for a given ver- 
tex and muon angle. The muon vertex, track angle, and 
energy, are found with a maximal likelihood fit, with the 
energy being determined from the total tank charge. The 
neutrino energy for both samples is reconstructed from 
the observed muon kinematics, treating the interaction as 
a 2-body collision and assuming that the target nucleon 
is at rest inside the nucleus: 

1 2mpE^, + ml- ml -ml 

^u^T^ , (1) 

^ TOp - + cosd^ JEf^ - ml 

Here is the mass of the proton, m^ is the mass of the 
muon, mi is the mass of the neutron in CCQE events 
and of the A(1232) in CCl7r+, 0^, is the reconstructed 
angle of the muon with respect to the beam axis (in the 
lab frame), and E^^ is the reconstructed muon energy. 

The distributions of signal events in neutrino energy 
are obtained through a two step process. First, the 
aforementioned cuts are applied to select the CCl7r+ and 
CCQE samples. These samples can be characterized by 
the cut efficiency (the fraction of signal events in the data 
set that pass the relevant cuts) and the signal fraction 
(the fraction of events in a given sample that are in fact 
signal events). Second, a Monte Carlo simulation (MC) 
is used to predict the signal fractions and cut efficiencies; 
these values are then used to correct the raw numbers of 
events passing cuts. 

For our primary measurement, we define CCItt"*" signal 
as events with exactly one /i~ and one tt"*" escaping the 
struck nucleus (which we call CCl7r+-like events) and 
CCQE signal as those with exactly one /x~ and no pions 
(CCQE- like events). Both event classes may include any 
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number of protons or neutrons, but no other types of 
hadrons. The observed cross section ratio is then defined 
as the ratio of CCl7r+-hke to CCQE-Hke events and thus 
has not been corrected for re-interactions in the struck 
nucleus. The signal fraction of the CClTr^-like (CCQE- 
Hke ) sample is predicted to be 92% (83%) and the cut 
efficiency is predicted to be 26% (38%) in 500 cm. Table |l] 
gives the composition of the CCl7r+-like and CCQE-like 
signal events in the MC. 



Process 


Fraction of CCln+- 
Hke events (%) 


Fraction of CCQE- 
hke events (%) 


CCItt"'' Resonant 


86.0 


9.4 


CCl7r+ Coherent 


6.3 


0.2 


CCQE 


2.4 


85.4 


Multi-pion 


2.5 


0.02 


CCItt" 


1.0 


2.5 


DIS 


0.2 


< 0.01 


Other 


1.6 


2.5 



TABLE I: Predicted composition of CClyr^-hke (one f^^ and 
one TT^ in the final state) and CCQE-like (one and no 
pions in the final state) events. 

To map reconstructed to true energy, we form a migra- 
tion matrix Aij representing the number of MC events in 
bin i of reconstructed energy and bin j of true energy. We 
then normalize each reconstructed energy bin to unity to 
obtain an unsmearing matrix. This is equivalent to a 
Bayesian approach discussed in [T^; it differs from the 
standard matrix inversion method in that the resulting 
unsmearing matrix is biased by the MC distribution used 
to generate it. We account for this in our uncertainties 
by including a variation in the MC distribution used to 
generate the matrix. Because we have good data/MC 
agreement, this effect is small. The advantage of this 
method is that it avoids the problems of numerical in- 
stability and the magnification of statistical errors which 
occur in matrix inversion. This unsmearing procedure 
also proved insensitive to variations in neutrino energy 
reconstruction, confirming that it performs as intended. 

With all the correction terms put together, the cross 
section ratio in each energy bin i is: 



■cuts,j 



(2) 



where subscript i runs over bins in true neutrino energy, 
subscript j indexes bins in reconstructed neutrino en- 
ergy, Nx-cuts denotes the number of events passing cuts 
for X = CCItt"*", CCQE, / denotes a signal fraction, e 
denotes a cut efficiency, and J7 is a neutrino energy un- 
smearing matrix that acts on a reconstructed distribution 
to return the true distribution. 

Figure [l] shows the observed CCl7r+-like to CCQE- 
like ratio extracted from the MiniBooNE data, including 
statistical and systematic uncertainties. 
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FIG. 1: Observed CCl7r"'"-like/CCQE-like cross section ratio 
on CH2, including both statistical and systematic uncertain- 
ties, compared with the MC prediction [6]. The data have 
not been corrected for hadronic re-interactions. 



The dominant systematic uncertainties on the cross 
section ratio arise from four sources: the neutrino flux 
(which largely cancels in the ratio), the neutrino inter- 
action cross sections (which affect the background pre- 
dictions), hadron re-interactions in the detector, and the 
detector simulation (which describes light propagation in 
the oil) . In the region of highest statistics (about 1 GeV) , 
there is roughly an 8% fractional error on the ratio result- 
ing from hadron re-scattering in the detector, 6% from 
neutrino cross sections, 4% from the detector simulation, 
2% from the neutrino flux, and 2% from the statistics of 
the two samples. 

In addition to these errors, an uncertainty on the 
dependence of the predicted CCl7r+ cross section is as- 
sessed based on comparison to MiniBooNE data. This 
contributes less than a 3% overall error to the measured 
ratio. Additional variations testing the sensitivity of the 
result to the event selection scheme, reconstruction algo- 
rithm, energy unsmearing method, and predicted 7r+ mo- 
mentum distribution in CCItt^ events are also included 
in the total uncertainty shown in Figure [T] Each of these 
added sources contributes a 1-2% uncertainty to the ratio 
in the region of highest statistics. 

Unlike the result presented in Figure [T] the ratio re- 
ported by all prior experimental measurements [2 IU |3] 
has been one in which the effects of final state inter- 
actions (FSI) in the target nucleus have been removed 
using MC. Solely for the purpose of comparison, we now 
extract a similarly corrected value. The FSI-corrected ra- 
tio is defined as the ratio of CCItt^ to CCQE events at 
the initial vertex and before any hadronic re-interactions. 
Thus, the signal fractions and cut efficiencies for the 
FSI-corrected ratio include corrections for intra-nuclear 
hadron re-scattering based on the MC's model for nu- 
clear effects. The measurement proceeds exactly as for 
the observed ratio (Figure [T|) , except that now we define 
CCl7r+ and CCQE, rather than CCl7r+-like and CCQE- 
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like, events as signal for the respective samples. With 
these definitions, the CCQE (CCl7r+ ) sample has a sig- 
nal fraction of 72% (87%) and a cut efficiency of 37% 
(20%) in 500 cm. The FSI-corrected ratio is shown in Fig- 
ure[2l The corrections for final state interactions have un- 
certainties associated with them, introducing additional 
systematic error to the cross section ratio. The fractional 
error on the ratio due to these corrections is roughly 6% 
in the region of highest statistics. 
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FIG. 2: FSI-corrected CCl7r+ to CCQE cross section ratio 
on CH2 compared with results from ANL (-D2) [I] and K2K 
(Cg-ffg) |3|. The data have been corrected for final state in- 
teractions and re-scaled for an isoscalar target. 

Here we limit our comparison to those experiments 
which reported both CCQE and CCl7r+ cross sections, 
using the same energy bins for each of these interac- 
tions, so as to facilitate comparison with our measured 
CCl7r+/CCQE ratio. Our result agrees with both ANL, 
which used a deuterium target, and K2K, which used 
CgHs (Fig. [2]). In order to make this comparison, the 
MiniBooNE and K2K results have been re-scaled to an 
isoscalar target. To perform this correction, we rescale 
the ratio by a factor of (1 — r)sp, where r is the ratio of 
neutrons to protons in the target and Sp is the fraction of 
7r+ production that is predicted (by MC) to occur on pro- 
tons. The resulting scaling factor is 0.80 for MiniBooNE; 
for K2K we use the factor of 0.89 provided in [3]. The 
results have not been corrected for their differing nuclear 
targets nor for the application of explicit invariant mass 
requirements (although the latter are similar). ANL used 
an explicit cut on invariant mass VF < 1.4 GeV. While 
no invariant mass cut is used in this analysis, the Mini- 
BooNE spectrum is such that CCItt"*" events occur only 
in the region W < 1.6 GeV; similarly, K2K's measure- 
ment covers the region W < 2 GeV 

The dominant reason for the difference between the 
ratios presented in Figures [l] and [2] is intra- nuclear pion 
absorption in CCItt"'" events, which cause these events to 
look CCQE-like. As a result of 7r+ absorption, a signifi- 
cant number of CCItt events appearing in the numera- 
tor in Figure [2] are in the denominator in Figure [T] Thus, 



the FSI-corrected ratio, shown in Figure [2I is 15% to 30% 
higher than the observed ratio in our energy range. 

In summary, MiniBooNE has measured the ratio of 
CCl7r+-like to CCQE-like events for neutrinos with en- 
ergy 0.4 GeV < < 2.4 GeV incident on CH2. This 
is the first time such a ratio has been reported. Addi- 
tionally, the ratio of the CCItt"'" and CCQE cross sec- 
tions at the initial vertex has been extracted using MC 
to remove the effects of final state interactions, in order 
to facilitate comparison with previous experimental mea- 
surements. The results are summarized in Table II. The 
measured ratios agree with prediction [6j [9] and previous 
data !U[5]. 



(GeV) 


CCItt+ZCCQE 
(FSI corrected) 


CClTr+-like/CCQE-hke 
(observed) 


0.45 ±0.05 


0.045 ±0.008 


0.036 ±0.005 


0.55 ±0.05 


0.130 ±0.018 


0.100 ±0.011 


0.65 ±0.05 


0.258 ±0.033 


0.191 ±0.019 


0.75 ±0.05 


0.381 ±0.047 


0.278 ±0.028 


0.85 ±0.05 


0.520 ±0.064 


0.371 ±0.040 


0.95 ±0.05 


0.656 ±0.082 


0.465 ±0.053 


1.05 ±0.05 


0.784 ±0.100 


0.551 ±0.066 


1.15 ±0.05 


0.855 ±0.114 


0.607 ±0.077 


1.25 ±0.05 


0.957 ±0.132 


0.677 ±0.091 


1.35 ±0.05 


0.985 ±0.141 


0.700 ±0.097 


1.5 ±0.1 


1.073 ±0.157 


0.777 ±0.109 


1.7 ±0.1 


1.233 ±0.207 


0.904 ±0.137 


2.1 ±0.3 


1.318 ±0.247 


1.022 ±0.161 



TABLE II: The MiniBooNE measured CClTr+ to CCQE (Fig- 
ure[2]| and CCItt "^'-like to CCQE-like (Figure[T]) cross section 
ratios on CH2 including all sources of statistical and system- 
atic uncertainty. 
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